Association between dietary acid load and metabolic health status in overweight and obese adolescents

The relationship between dietary acid load (DAL) and metabolic health status in adolescents has not been studied yet. We aimed to examine the association between DAL and metabolic health status in Iranian overweight/obese adolescents. This cross-sectional study included 203 overweight/obese adolescents selected by a multistage cluster random sampling method. Dietary intakes were assessed using a validated 147-item food frequency questionnaire (FFQ). Anthropometric indices and blood pressure values were measured. Fasting blood samples were obtained to determine glucose, insulin, and lipid profiles. Based on two methods (International Diabetes Federation (IDF) criteria and combination of IDF with Homeostasis Model Assessment Insulin Resistance (HOMA-IR)), participants were classified into metabolically healthy obese (MHO) or unhealthy obese (MUO). Adolescents in the highest tertile of potential renal acid load (PRAL) and net endogenous acid production (NEAP), compared with those in the lowest tertile, had 172% (95% CI 1.32–5.59) and 161% (95% CI 1.26–5.41) higher odds of MUO status, based on IDF criteria. This association was significant after adjustment for age, sex, and energy intake (PRAL: OR 2.42; 95% CI CI 1.13–5.15; NEAP: OR 2.52; 95% CI 1.17–5.41); but it disappeared after adjustment for other confounders. Based on IDF/HOMA-IR definition, there was a significant positive association between PRAL and being MUO only in the crude model (OR 2.37; 95% CI 1.13–4.96). The stratified analysis revealed that these associations for NEAP scores were stronger among overweight subjects than obese individuals, based on both metabolic status definitions. However, after adjustment for all potential confounders these relations were insignificant. Having higher DAL might be associated with higher odds of MUO phenotype in Iranian overweight/obese adolescents. More prospective studies are warranted to confirm this finding.


Methods
Study design and population. This cross-sectional study included 102 girls and 101 boys (12-< 18 years old) who were selected by a multistage cluster random sampling method. Participants were randomly selected from 6 different districts of the city of Isfahan, Iran, using a stratified, multi-stage cluster sampling design. In total, 16 schools were chosen. Body weight and height of all students were measured and BMI was calculated based on weight (kg)/height 2 (m). Then, adolescents were classified into 3 groups: normal-weight, overweight, and obese 25 . Overweight and obese students with different socioeconomic status (SES) were considered to participate in the study. The students who had the following conditions were not included in this study: having genetic or endocrine disorders (such as type 1 diabetes mellitus, hypothyroidism, and Cushing's syndrome), being on a weight-loss diet, taking vitamin and mineral supplements or medications that might affect body weight, hypertension, blood glucose or lipid profile. All participants and their parents signed the written informed consent. The study protocol was approved by the local Ethics Committee of Isfahan University of Medical Sciences. www.nature.com/scientificreports/ Assessment of dietary intakes. Dietary intake was assessed using a validated 147-item semi-quantitative food frequency questionnaire (FFQ) 26 . According to this questionnaire, dietary intake in the previous year was calculated by determining the frequency of consumption of each food item on a daily, weekly, or monthly basis. The amount of food consumed was evaluated by asking a standard common portion size. We converted the portion sizes of consumed foods to grams per day using household measures. Finally, the grams of food intakes were entered into Nutritionist IV software to derive nutrient intake data. Nutritionist IV software was based on USDA food composition data; some Iranian foods were also incorporated into it.
Calculation of DAL indices. Urinary net acid excretion depends on dietary intake, and direct measurement is difficult. Therefore, two dietary indices (PRAL and NEAP) have been defined to calculate body acid production or DAL. In a previous study by Remer et al., the validity of both PRAL and NEAP scores compared to a 24-h urinary acid load was confirmed 27 . The PRAL score was calculated based on the following formula 27 : In addition, we estimated the NEAP score as follows 28 :

Assessment of anthropometric indices and cardiometabolic risk factors. Two trained nutrition-
ists measured weight using a calibrated electronic scale (Seca Instruments, Germany) in minimal clothing and without shoes to the nearest 0.1 kg. Also, standing height was measured without shoes with a stadiometer to the nearest 0.1 cm. Then, body mass index was calculated as kg/m 2 . Based on the sex-and age-specific World Health Organization (WHO) cut-off points, participants were classified as overweight or obese adolescents 25 . Waist circumference was recorded twice at the midway between the lowest rib and the superior border of the iliac crest by the use of an unstretched flexible anthropometric tape at 0.1 cm precision, after a normal expiration and without any pressure on the body surface. The mean of the two measurements for each student was considered as WC. Systolic blood pressure (SBP) and diastolic blood pressure (DBP) were measured twice at the right arm using a mercury sphygmomanometer with a suitable cuff size after a rest period of 15 min. We considered the average of the two measurements in the analysis. Venous blood samples after 12 h of fasting were collected to determine biochemical values. Fasting blood glucose (FBG), high-density lipoprotein cholesterol (HDL-c), triglycerides (TG), and insulin concentrations were measured with standard methods. Homeostasis Model Assessment Insulin Resistance (HOMA-IR) was calculated to estimate Insulin resistance (IR) by using the following formula: HOMA-IR = (insulin (µUI/mL) × glucose (mg/dL))/405. Assessment of metabolic status. Two methods were used to classify participants into MHO or MUO individuals. Students with two or more of the four following risk factors were defined as MUO subjects, based on the modified International Diabetes Federation (IDF) criteria 29 : (1) increased fasting blood glucose (≥ 100 mg/ dL), (2) increased blood pressure (≥ 130/85 mmHg), (3) increased triglycerides (≥ 150 mg/dL) and (4) decreased HDL-c (< 40 mg/dL for the age of < 16 y, and < 50 mg/dL in girls/ < 40 mg/dL in boys for the age of ≥ 16 y). MHO adolescents were those with one or no risk factors. In the second method, presence of insulin resistance (defined by HOMA-IR score) was added to the IDF criteria that were used in the first definition 30 . Students with HOMA-IR score ≥ 3.16 that had two or more above-mentioned metabolic risk factors were considered as MUO individuals and those with HOMA-IR < 3.16 were classified considered as MHO adolescents. Based on some previous studies on children and adolescents, the cut-off value for HOMA-IR was set at 3.16 31-33 . Assessment of other variables. The Physical Activity Questionnaire for Adolescents (PAQ-A), which contains 9 questions on various aspects of physical activity over the past week, was used to determine the physical activity level of participants 34 . Items 1 to 8 of this questionnaire were scored from 1 to 5, with 1 representing the least physical activity and 5 representing the most. The ninth question examined the unusual activity of adolescents during the previous 7 days. Based on their total scores, adolescents were categorized as highly active (score ≥ 4), moderately active (4 < score ≤ 3), lowly active (3 < score ≤ 2), sedentary, or not having an orderly week of physical activity (score < 2). A validated demographic questionnaire was used by trained investigators to ascertain SES of students (29), based on the following items: parental education level, parental job, number of family members, having cars in the family, having personal room, having computers/laptops and taking trips in the year. Finally, a total score of SES (ranged between 4 and 27) was calculated. Information of sex, age, history of diseases, and use of medications and dietary supplements of participants was gathered through a questionnaire.
Statistical analysis. We presented continuous variables as mean ± SD/SE and qualitative variables as frequency (percentage). Participants were categorized into tertiles of PRAL and NEAP (based on their scores). We used one-way analysis of variance (ANOVA) and the χ2 test to compare quantitative and categorical variables between tertiles of PRAL and NEAP. Age-, sex-and energy-adjusted dietary macro and micronutrient intakes of participants were compared across tertiles of PRAL and NEAP using analysis of covariance (ANCOVA Ethical approval and consent to participate. The study procedure was performed according to declaration of Helsinki and STROBE checklist. All participants provided informed written consent. The study protocol was approved by the local Ethics Committee of Isfahan University of Medical Sciences.

Results
The present study included 101 (49.8%) boys and 102 (50.2%) girls with a mean age of 13.34 ± 1.38 and 14.61 ± 1.57 and BMI of 26.67 ± 2.82 and 28.03 ± 3.48, respectively. General characteristics and cardiometabolic risk factors of study participants across tertiles of PRAL and NEAP are presented in Table 1. Adolescents in the highest tertile of PRAL in comparison to those in the lowest tertile had higher FBG and TG (P < 0.05). Furthermore, subjects in the top category of PRAL had lower physical activity level and different SES status (P < 0.05). Adolescents in the highest tertile of NEAP had higher FBG than those in the lowest tertile and were different in terms of physical activity from those in the first category (P < 0.05). There were no significant differences in other general characteristics and cardiometabolic risk factors among tertiles of PRAL and NEAP. Dietary intakes of participants across tertiles of PRAL and NEAP are presented in Table 2. Participants in the top tertile of PRAL, compared to the bottom tertile, had greater intake of energy, protein, thiamin, niacin, selenium and sodium and lower intake of fat, MUFA, vitamin C, vitamin A, riboflavin, vitamin B 6 , folate, vitamin Table 1. General characteristics and cardiometabolic factors of study participants across tertiles of PRAL and NEAP 1 . 1 Values are Mean ± SD; unless indicated. Abbreviations: NEAP: net endogenous acid production; PRAL: potential renal acid load; BMI: Body Mass Index; HOMA-IR: Homeostasis Model Assessment Insulin Resistance; HDL-c: high-density lipoprotein cholesterol; 2 P value for one-way analysis of variance (ANOVA) and χ2 test for quantitative and categorical variables, respectively. 3 Socioeconomic status (SES) score was evaluated based on parental education level, parental job, number of family members, having car in the family, having computer/laptop, having personal room and having travel by using demographic questionnaire.

Discussion
We found significant direct associations between PRAL and NEAP with odds of MUO status in Iranian adolescents, based on IDF criteria. However, these associations were dependant to covariates. The same finding was observed between PRAL and odds of MUO based on IDF/HOMA-IR definition, before adjustments for potential confounders. Stratified analysis indicated that greater NEAP categories were associated with higher odds for MUO among overweight adolescents, based on both definitions of metabolic health status. The present study is one of the first investigations on the association between DAL and metabolic health status in adolescents. The attainment of MHO may be a worthwhile first goal in the management of obesity and its complications, although MHO is an intermediate state before the development of MUO. Diet is one of the most important factors in determining obesity phenotype that should be considered to improve general health status 13 .
Our study showed possible beneficial relations between lower acid load of diet and metabolic health. Some prior studies have also investigated the association between DAL and obesity or metabolic syndrome. It is also worth noting that although some components of MHO/MUO and metabolic syndrome were the same, these are two different outcomes. Arisawa et al. have reported significant positive trends between NEAP scores with metabolic syndrome and some components of MetS in adults 35 . A systematic review and meta-analysis on observational studies additionally revealed that high DAL content was associated with higher obesity prevalence and serum TG concentrations 17 . A recent study demonstrated that the odds of obesity and abdominal adiposity significantly increased in tertile of DAL and also DAL could be an indicator of diet quality 36 . However, a previous cross-sectional study by Rahbarinejad et al. indicated no significant association between PRAL and NEAP and MetS in children before or after adjustment for potential confounders 24 . Also, another study found the same finding in adults, but higher scores of NEAP were associated with increased risk of impaired fasting blood glucose 37 . A national multi-center cross-sectional survey on 5326 Iranian students aged 6-18 years after evaluating anthropometric indices, showed a direct association between diet-induced acid load with neck circumference 23 . Furthermore, an inverse relationship was found between DAL indices and parental BMI 23 . www.nature.com/scientificreports/ www.nature.com/scientificreports/ We observed that greater scores of PRAL and NEAP were related to higher FBG. Also, individuals with higher score of PRAL had higher TG. Some other investigations have also evaluated the relationship between DAL and cardiometabolic risk factors. A cross-sectional study on 371 Iranian adult women pointed out that women with higher NEAP scores had higher weight, WC, and TG concentrations. In case of PRAL, only a significant relation with TG was found 20 , which was similar to our finding. Lee et al. reported that diets with high acid-forming potential (high PRAL and NEAP scores) were positively associated with higher risk of IR during 7.3 years followup in a large population-based cohort study 38 . In another study on Japanese workers, high DAL was associated with IR, but not with FBG or HbA1c levels 21 . Another study by Haghighatdoost et al. demonstrated that PRAL and protein to potassium ratio were directly related to HbA1c among diabetic nephropathy patients 39 . Dehghan et al. revealed a potent negative relation between high DAL, particularly with PRAL scores, and cardiometabolic risk factors (including SBP, DBP, insulin, and odds of diabetes) in a systematic review and meta-analysis 40 . Differences in the findings of previous studies might be associated to the study populations, dietary intake assessment tools, sample size of the studies, applied methods to calculate dietary acid load, definition of dietary patterns, and various controlled confounders.
Several possible mechanisms might account for the relationship between DAL and MUO. Our study showed that adolescents in the top tertiles of PRAL and NEAP had greater intake of protein and carbohydrate, respectively. Animal food sources such as meat, fish, egg, chicken, cheese, and also cereals are rich in sulfur-containing amino acids, phosphorous and chloride. Consumption of these foods can contribute to higher acid production in the body; while high content of malate, citrate and glutamate in fruits and vegetables can contribute to base production. Therefore, animal-based foods, which are high in western diets, can be considered the most acidproducing diets and are associated with metabolic diseases 41 . These foods are pro-inflammatory and can increase inflammatory cytokines 42 , which are involved in lipid disorders and IR, particularly in obese subjects 43 . Higher DAL stimulates the production of glucocorticoids 44 . High cortisol levels induce lipase activity (lipoprotein lipase and hormone-sensitive lipase), which can lead to the increased efflux of fatty acids to the bloodstream and augment the production of very-low-density lipoproteins in the liver 45 . Furthermore, increased acidity of the body is related to the renal excretion of magnesium, calcium, and potassium 46 , although we did not measure serum concentration of these nutrients in our study. Increased BP and IR and obesity may result from these processes 47,48 . Also, diet-induced acidosis causes a decrease in adipokines (leptin and adiponectin), which can impair appetite suppression 49,50 . In addition, it has been demonstrated that high DAL reduces lean body mass and ultimately leads to increased body fat synthesis 51 .
The present study has some strength. First, the novelty of this study was the evaluation of DAL-metabolic health relation in a sample of adolescents from all socioeconomic districts of a large central city in Iran. Second, two different definitions of metabolic health status were used to determine MUO/MHO status. Third, several potential confounders were adjusted in the analysis. However, we acknowledge some limitations. The cross-sectional design of the study did not allow us to examine causality; so, conducting prospective studies is necessary. In addition, recall bias and other potential reporting biases may have affected the results. Also, the interview setting was used to collect dietary intakes data, which may be subject to social desirability biases. Despite controlling many potential confounders, residual confounders (such as the pattern of dietary habits, puberty, and sleep health) may affect the association between DAL scores and MUO. Finally, although BMI and WC were measured to define obesity and abdominal obesity, we could not measure body composition and fat distribution that could be involved in metabolic health status.
In conclusion, our findings indicated that diet-induced acid load might be associated with higher odds of MUO phenotype in overweight/obese adolescents. Thus, lower acidogenic food ingredients in the diets and high consumption of fruits and vegetables may improve the metabolic health status of adolescents. Further studies, particularly prospective nature, are needed to confirm these findings.

Data availability
The datasets generated and/or analysed during the current study are not publicly available as per the rules and regulations of the Isfahan University of Medical Science but are available upon reasonable request from the corresponding author. www.nature.com/scientificreports/ 51. Faure, A. M., Fischer, K., Dawson-Hughes, B., Egli, A. & Bischoff-Ferrari, H. A. Gender-specific association between dietary acid load and total lean body mass and its dependency on protein intake in seniors. Osteoporos. Int. 28, 3451-3462. https:// doi. org/ 10. 1007/ s00198-017-4220-z (2017).